Superfluid pairing between fermions with unequal masses 
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We consider a superfluid state in a two-component gas of fermionic atoms with equal densities 
and unequal masses in the BCS limit. We develop a perturbation theory along the lines proposed 
by Gorkov and Melik-Barkhudarov and find that for a large difference in the masses of heavy 
(M) and light (m) atoms one has to take into account both the second-order and third-order 
contributions. The result for the critical temperature and order parameter is then quite different 
Ch ' from the prediction of the simple BCS approach. Moreover, the small parameter of the theory 

^ I turns out to be (j)F\a\/h)y^ M/m <^ 1, where pp is the Fermi momentum, and a the scattering 

length. Thus, for a large mass ratio M/m the conventional perturbation theory requires significantly 
smaller Fermi momenta (densities) or scattering lengths than in the case of M ~ m, where the small 
parameter is pF\a\/h <^ 1. We show that 3-body scattering resonances appearing at a large mass 
ratio due to the presence of 3-body bound Efimov states do not influence the result, which in this 
, ^ ■ sense becomes universal. 

D 

r-| PACS numbers: 03.75.Ss, 74.20.Fg 



I. INTRODUCTION 



Superfluid pairing in a two-component gas of fermions is a well-known problem [l| lying in the background of 

■ extensive studies in condensed matter and nuclear physics 0, H, 3 ■ Recently, this problem was actively investigated 
Ch ' in cold gases of fermionic atoms (see [H for a review). Experimental efforts were focused on ^Li or *'^K atoms in two 
^ • different internal (hyperfine) states, where one can use Feshbach resonances for switching the sign and tuning the 

absolute value of the interspecies interaction (scattering length a), which at resonance changes from — oo to -foo. In 
this respect, one encounters the problem of BCS-BEC crossover discussed earlier in the context of superconductivity 
' BBS [13 and for superfluidity of two-dimensional ■'He films On the negative side of the resonance (a < 0), 

one should have the Bardeen-Cooper-Schrieffer (BCS) superfluid pairing at sufficiently low temperatures, and on the 
positive side (a > 0) one expects Bose-Einstein condensation of diatomic molecules formed by atoms of different 
QQ , components. Remarkable achievements of cold-atom physics in the last years include the observation of superfiuid 
behavior through vortex formation in the strongly interacting regime (n |a|'^ > 1, where n is the gas density) |13| . and 
the formation and Bose-Einstein condensation of long-lived weakly bound diatomic molecules at a > [iJI . Ongoing 
experiments with atomic Fermi gases have reached temperatures in the nanokelvin regime, where at achieved densities 

■ one has T ^ O.lEp, with Ep being the Fermi energy. For a < the experiments are now approaching superfluidity 
, in the BCS limit where n\a\^ ^ 1. 

IL^ ■ Currently, a new generation of experiments is being set up. In particular, it is dealing with mixtures of different 
fermionic atoms or mixtures of fermions and bosons. The main goal is to reveal the infiuence of the mass difference 
on superfluid properties and to search for novel types of superfluid pairing. The first experiments demonstrating a 
possibility of using Feshbach resonances and creating coUisionally stable mixtures of with ^Li and/or with *^Rb, 



9^ 1 and ^Li with ^^Na have already been performed [ll, [T^, [l3, [3, E^, [20II . Recent theoretical literature on mixtures 
of different fermionic atoms contains a discussion of the BCS limit l2ll . [2^ . [23| , the limit of molecular BEC [IJ , 
BCS-BEC crossover [l^H^H^I, and the strongly interacting regime j28| . 

In this paper we consider a two-component mixture of fermionic atoms with different masses and attractive inter- 
component interaction in the BCS limit. It is assumed that the densities of the two species are equal which means that 
there is no mismatch between their Fermi surfaces, leading to the usual BCS type of superfluid pairing. Other kinds 
of pairing that can occur and compete with BCS, especially for unequal densities, will be discussed elsewhere [29j . 
Here, we generalize the perturbation treatment of the gap equation, introduced by Gorkov and Melik-Barkhudarov 
[sot for equal masses of fermions belonging to different components. This approach takes into account the interaction 
between the atoms in a Cooper pair due to the polarization of the medium and allows one to correctly determine the 
dependence of the zero-temperature gap Aq and superfluid transition temperature Tc on the masses of heavy (M) 
and light (m) fermionic atoms. As we shall see below, already the second order of the perturbation, the so-called 
Gorkov-Mclik-Barkhudarov contribution, leads to a very different dependence of the preexponcntial factor in the 
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expressions for Aq and Tc on the mass ratio Mjm. compared to the prediction of the simple BCS theory. 

For a large mass ratio Mjm ^ 1, we include higher order contributions and show that the actual small parameter 
of the theory is {pf\o\IK) ^Mjm ^ 1 (pi? is the Fermi momentum), not simply pi?|a|/fi, <C 1 as in the case of M ^ to. 
We give a physical interpretation of this fact and calculate effective masses of heavy and light fermions. 

Large mass ratios Mjm are realized in electron-ion plasmas, where the heavy ion component is usually considered 
as non-degenerate [3l| . The electron-proton pairing in the hydrogen plasma, assuming quantum degeneracy for both 
electrons and protons, was discussed by Moulopoulos and Ashcroft |32l |. They found that at low temperatures the 
Coulomb electron-proton attraction leads to the appearance of a (momentum-dependent) gap which for sufficiently 
high densities is comparable with the Coulomb interaction at the mean interparticle separation. Note that this problem 
is quite different from ours where the attractive interaction between heavy and light fermions is short-ranged. 

Before proceeding with our analysis we make two important remarks. First of all, if the masses of heavy and light 
fermionic atoms are very different from each other and the mass ratio exceeds a critical value, M/m > 13.6, then 
two heavy and one light fermion can form 3-body weakly bound states. The appearance of these states, predicted 
by Efimov [33| . can be easily understood in the Born-Oppenheimer picture (40| . If we fix the two heavy atoms at a 
relative distance R < |a|, a localized state for the light atom appears due to the presence of the heavy pair, which in 
turn mediates an attractive interaction ~ —h^/mR^ between the heavy atoms (see, e.g. [l^l and references therein). 
For a large mass ratio, M/m > 13.6, this mediated attraction overcomes the kinetic energy of the relative motion of 
the heavy atoms and one has the well-known phenomenon of "fall into center" [IJl . The energy of this state is bounded 
from below only due to short-range repulsion. The corresponding wave function of the relative motion of heavy atoms 
acquires a large number of nodes thus showing the presence of many bound states. This makes the 3-body problem 
non-universal in the sense that aside from the 2-body scattering length a, the description of this problem requires 
one more parameter - the so-called 3-body parameter coming from short-range physics. Also, the presence of weakly 
bound Efimov states introduces a resonant character to the 3-body scattering problem. This is especially important 
for the Gorkov-Melikh-Barkhudarov contribution as it is actually dealing with processes involving 3 particles. We, 
however, have found that the 3-body resonances are rather narrow and their contribution is not important. This 
makes the Gorkov-Melikh-Barkhudarov approach universal at any mass ratio M/m. 

Our second remark is related to analogies between BCS pairing of unequal-mass particles in cold-atom and high 
energy physics. We wish to emphasize that there are strong physical differences between the pairing of particles of 
different masses in relativistic and in nonrelativistic systems such as cold atoms. The problem arises in relativistic 
systems in the study of hadronic matter (35| . It is thought that, at the high densities achieved in neutron stars, quarks 
become deconfined and the different types of them (e.g. up, down and strange quarks) will tend to form Cooper pairs 
with each other. These different types of quarks have different masses. This relativistic limit has been investigated by 
Kundu and Rajagopal [36| and is characterized by the Fermi momentum being larger than the bare mass: pp 3> mc. 
Linearizing the momentum near the Fermi surface p = pp + Sp, {Sp/pp <^ 1), we can expand the free particle energy 
to the lowest nonvanishing order in Sp/pp and mc/pp: 



We see that, as far as the kinetic energy is concerned, a change in mass will amount to a shift in the chemical potential 
of the species which is proportional to to^ and depends inversely on pp. Therefore, pairing between particles with 
different masses in the relativistic limit is equivalent to studying the problem of pairing of equal mass particles in the 
presence of a difference between the chemical potentials of the two species. In the nonrelativistic limit (pp <^ mc) 
the situation is different: 



and so the mass change cannot be incorporated into the chemical potential, requiring a very different analysis which 
we carry out here. 

The paper is organized as follows. In Section|TT]we present general equations, and in Sections III and IV we calculate 
the critical temperature Tc. Section V is dedicated to the discussion of the small parameter of the theory, and in 
Section VI we discuss the order parameter and excitation spectrum. In Sec. VII we analyze the three-body resonances 
at a large mass ratio M/m and show that they do not change the result of the Gorkov-Melik-Barkhudarov approach. 
In Sec. VIII we conclude. 
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V^>(P+P') = 




Figure 1: The leading many-body contribution (second order, or Gorkov-Melik-Barkhudarov correction) to the effective in- 
teraction between heavy (thick line) and light (thin line) fermions. The dashed line corresponds to the coupling constant 



II. GENERAL EQUATIONS 



We consider a uniform gas composed of heavy and light fermionic atoms with masses M and m, respectively. Both 
heavy and light atoms are in a single hyperfine state, and considering low temperatures we omit heavy-heavy and 
light-light interactions. The interaction of heavy atoms with light ones is assumed to be attractive and characterized 
by a negative s-wave scattering length a < 0. The Hamiltonian of the system has the form: 



H ^ dr 



i=1.2 ^ ' 



A (r) + g-ipt i^)i^2 (r)^2 (r) 



(1) 



where 4>i{r) and tp^' arc the field operators of fermionic atoms labeled by indices i — 1 (heavy) and « = 2 (light), jii 
is the corresponding chemical potential, and g = 27rh?a/mr is the coupling constant, with nir — Mm/ {M + m) being 
the reduced mass. Since the densities of the two species are equal, ni = n2 = n, they have the same Fermi momentum 
Pf ~ h{6TT^ny^^ , and, hence, = p\/2M and ^2 = Pp/2m. Finally, we require that the system be in the weakly 
interacting regime, which requires the inequality pp\a\/h <^ 1. 

We now consider the usual BCS scheme where a heavy atom with momentum p is paired to a light one having 
momentum —p. This leads to the gap equation 



,M-/(i?+(p'))-/(i?_(p')) 



E+{p') + E-ip') 



A(p') 



(2) 



where f{E) = [exp(E /T) + 1] ^ is the Fermi-Dirac distribution function, and we assume that the order parameter is 
real. The dispersion relations for the two branches of single-particle excitations are written as 



E±ip) = ± 



6(p)-e2(p 




6(p)+e2(p) 



(3) 



and the quantities ^1,2 are given by ^i(p) = {p^ —pp)/2M and ^2(p) = {p^ — pp)/2m. The function yoff(p,p') = 
g + 6V{p, p') is an effective interaction between particles in the medium, where the quantity 6V{p, p') originates from 
many-body effects and is a correction to the bare interparticle interaction g. The leading correction is second order 
in g and the corresponding diagram is shown in Fig. [T] 

The integral in Eq. ([2|) diverges at large momenta due to the first term in T4ff . This divergency can be eliminated 
by expressing the bare interaction g in terms of the scattering length a P, Isol . [3^ . If we confine ourselves to the 
second order in perturbation theory with respect to g, then the renormalized gap equation reads 



A(p) 



dp' 

{2Trhy 



dp' 

{2TThf 



■5V{p,p' 



l-f{E+{p'))-f{E^{p')) 2mr 



E+{p')- 
fiE+{p')) 



E-iP') 
-f{E-{p')) 



A(p') 



E+{p') + E^ip') 



A(p') 



(4) 
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The convergence of the integral over p' in the first term of the right-hand side of Eq. (|4]) is now obvious, while the 
convergence of the second term is due to the decay of SV at large momenta (sec Eq. pT|) below) . The gap equation 
accounting for higher orders in g will be derived and discussed in Section IV. 
In the limit of A — > one can reduce Eq. ([¥]) to a linearized gap equation: 



A(p) = - 



2Trh'^a 
dp' 

{2'Khy 



dp' 

{2-Khf 



tanh[a(p')/2r] + tanh[e2(p')/2r] 2? 



2Ki(p')+6(p')] 



A(p') 



,,r, „,, tanh[ei(p0/2r] + tanhfc(pO/2T] ^ 

^'^^^P' P ^ 2[e,(p') + 6(P')] ^ 



(5) 



The critical temperature Tc is determined from Eq. ([5]) as the highest temperature at which this equation has a 
non-trivial solution for A. 



III. CRITICAL TEMPERATURE. BCS AND GM APPROACHES 

The first line of Eq. (O corresponds to the linearized gap equation in the traditional BCS approach: 



A(p) = - 



dp' 



(27rfi)3 



tanh[^i(p')/2T] + tanh[^2(p')/2r] 2?- 



2Ki(p') + 6(p') 



A(p') 



(6) 



In this case the order parameter is momentum independent, A(p) = A, and Eq. ^ reduces to the equation for the 
critical temperature: 



1^^ 
2 



In 



8/ii exp(7 - 2) 



In 



8/^2 exp(7 - 2) 



ttTbcs t^Tbcs 
where 7 = 0.5772 is the Euler constant and we introduced a small parameter 

A = 2pF \a\ /ttH < 1. 

Equation ([7]) is obtained straightforwardly. First, integrating Eq. ([5]) over the angles one has 



(7) 



(8) 



1 = - / dx [tanh[(a;2 - 1)mi/2T;] + tSinh[{x^ - l)/i2/2Tj - 2] 
2 Jo 

tanh[(x2 - l)^i/2rj + tanh[(a;2 - 1)^2/2TJ 



dx- 



wherc x = p/pf- For fii/Tc >> 1 the integrand of the first integral is equal to —4 for a; < 1, and for a; > 1 it rapidly 
drops to in a narrow interval of x, where |a; — 1| < Tc/ fJ. << 1. The contribution of this interval can be neglected 
and, therefore, the first term equals — 2A. Then, after integrating the second term by parts, we obtain 



dx 



C/U2/T'c 



cosh2[(a;2 - l)/ii/2rj cosh2[(x2 - l)/i2/2rj_ 



In 


x-l 






} 




x + 1 





The final integration can easily be performed by using the fact that the integrand is non-zero only in a narrow range 
of X, where |x — 1| < Tc/fi << 1. We can therefore introduce a new variable y ~ x — I and extend the limits of 
integration over y from —00 to -l-oo. The equation then reads: 



1 = A 



dy 



cosh^(j//ii/rc) cosh^(j//i2/Tc) 



In- 



and performing the integration one arrives at Eq. ([T]). This equation gives the critical BCS temperature (cf. [22 
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Tbcs = - exp(7 - 2)y/ij.ifj,2 cxp y-jj ■ (9) 

However, the linearized BCS gap equation ^ can only be used for the calculation of the leading contribution to 
the critical temperature, corresponding to the term ~ in the exponent of Eq. ([9]). Therefore, only the exponent 
in this equation is correct. As was shown by Gorkov and Melik-Barkhudarov [lOl, the preexponential factor in Eq. ^ 
is determined by next-to-leading order terms, which depend on many-body effects in the interparticle interaction. 
These are the interactions between particles in a many-body system through the polarization of the medium - virtual 
creation of particle-hole pairs. 

The importance of the many-body effects for the preexponential factor can be understood as follows. After per- 
forming the integration over momenta, the gap equation ([5]) can be qualitatively written as 



1 = i^fV,s ^n^ + C , (10) 

where = rrirPF /t^'^^i^ ■ In this formula, the large logarithm In/^/Tc comes from the integration over momenta near 
the Fermi surface, whereas the momenta far from the Fermi surface contribute to the constant C which is of the order 
of unity. We then write lypVcs = — A + aA^, where the first term is the direct interparticle interaction and we keep only 
the second order term in the many-body part SV of the effective interaction. It is now easy to see that Aln/i/Tc 1 
and, therefore, the terms aX^ In fx/Tc and AC are of the same order of magnitude. As a result, both terms have to be 
taken into account for the calculation of the preexponential factor. Also, note that the contribution of the many-body 
part of the interparticle interaction comes from momenta near the Fermi surface, which are responsible for the large 
logarithm In^/Tc ~ A""'^. Hence, only the values of SV at the Fermi surface are important. 

We now calculate the contribution of the many-body effects to the preexponential factor for the critical temperature. 
They are usually called Gorkov-Melik-Barkhudarov (CM) corrections. As it was argued above, in the weak coupling 
limit the most important contributions to the effective interaction are second order in g (the role of high order terms 
will be discussed later). In the considered case of a two-component Fermi gas with an s-wavc interaction, there is 
only the contribution shown in Fig. [U and the corresponding analytical expression reads: 



6V{p,p') 



dk /[ei(k + q/2)]-/K2(k-q/2)] 
(27r)3 6(k + q/2)-6(k-q/2) ' 



(11) 



where q = p + p'. In obtaining this expression we used the zero-temperature distribution function /[$i,2(p)] ~ 
0(— ^1 2(p)), with 9{x) being the step function. This is legitimate because the finite temperature corrections are 
proportional to the ratio of the critical temperature to the chemical potential and, therefore, are exponentially small. 
As can be seen from Eq. pT|) . the effective interaction SV{p, p') changes on the momentum scale p ^ p' ^ pp. 

We now solve Eq. ([5]). In this equation, the momentum dependence of the order parameter originates only from 
the momentum dependence of the many-body contribution to the interparticle interaction and, therefore, contains 
an extra power of the small parameter A. As a result, this dependence can be ignored in the first integral on the 
right-hand side of Eq. (O, and we can simply replace there the order parameter A(p') by its value on the Fermi 
surface A{pf). This does not affect the convergence of the integral at large momenta and, hence, changes only the 
constant C in Eq. pop . The corresponding modification, however, is proportional to the small parameter A and can 
be neglected. In the second integral on the right-hand side of Eq. ([S]), as we have discussed earlier, only momenta 
p' near the Fermi surface {p' « pp) are important, and we can also put p' = pp in A{p') and SV{p,p'). The gap 
equation then reads: 



A(p) 



27rft2a 



dp' 
(2^ 



tanh[Ci(p')/2rc] + tanh[C2(p')/2rc] 2r7v 



p<ApF 



dp' 



■5V{p, n'pp 



2[a(p')+6(p')] / 

tanh[ei(p0/2r,] + tanh[e2(p0/2re 

2Ki(P') + e2(p')] 



A{n'pp) 



A(nV), 



(12) 



where n' is the unit vector in the direction of p' and we introduced an upper cut-off App with A 1, for the 
purpose of convergence at large momenta. The exact value of A is not important because, as we mentioned above, the 
contribution of large momenta to this integral has to be neglected. To derive an equation for the critical temperature, 
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we consider Eq. ([T2|) for p = npp and average it over the directions of n. Taking into account that the order parameter 
for the s-wave pairing can only depend on the absohite value of the momentum, we obtain: 



1 = 



where 



dp' 

(2^ 



tanh[^i (p') /2TJ + tanh[6 {p')/2Tc] 2m,. 



2Ki(p')+6(p')] 



n'2 



^6V 
= Aln 



dp' tanh[a(p')/2Td + tanh[e2(p')/2Tj 



8^^ijU2exp(7 - 2) 



2Ki(p') + 6(p')] 

jTj, 87A'i7*2exp(7 + A-2) 



— vp5V In 



(13) 



(14) 



is the s-wave component of the many-body interaction. Using Eq. (fTTj) and integrating over the angles in Eq. (|14p . 
we obtain 



(15) 



where k = Mjm and the function /(k) is given by 



/(k)^- 3(1 + At) /■'^ 1^ 

^ 4(l + ln4)io 7o 

A straightforward lengthy integration of Eq. (fT6|) yields 



(p2-l)(K-l) + 4g(p-g) 



(p2-l)(K-l)-4q(p + (?) 



(16) 



(1 + k) 



4(1 + In 4) 
4 / (k 3) 



K + 1 



3k 

K + 2 



ln(2) 



K - 1 



In- 



3k 

K + 1 



■Ink- 11 



3(k-1) V6(k-1)2 6^ y 2 
From Eqs. (|13p and (jl5p we obtain the following expression for the critical temperature: 



(17) 



Tgm = -eT-27^exp[-(A - ^.^JV^)-!] « Tbcs^M~vfW 

TT 



eT /'2y^^^^p| l + ln4 



[/(k) + /(k-1) - 1]| V^exp (^-i^ 



It can be rewritten in the form 



^GM = ^ f-j Mi^^(«)cxp f--j = 0.277 ^F(K)exp(-7rft/2pi.|a| 



(18) 



(19) 



where we expressed the Fermi energy of light atoms [I2 through the heavy-atom Fermi energy /ii ~ p^p/2M, and the 
function F{k) is given by 



F{k) ~ \/kcxp < — 



1 + ln4 



[/(k) + /(k-1)-1]| 



(20) 



For equal masses one has k = M /m = 1, and Eqs. p7p and (PD|) give /(I) = 1/2, F{1) ~ 1. Then Eq. (fT5|) reproduces 
the original result of Ref. [s^l. The function F(k) is shown in Fig. [2l For a large mass ratio k ^ M/m ^ I, this 
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Figure 2: The function F{k) in the preexponential factor of Eq. (|19|l . 





«Li 


«K 


»^Sr 


"^Yb 


«Li 


1.000 


0300 


0.161 


0.090 


«K 




0.150 


0.097 


0.062 


»^Sr 






0.069 


0.048 


Yb 








0.035 



Table I; Critical temperatures Tqai for mixtures of various atomic species as given by Eq. (|19|) . The scattering length a and 
Fermi momentum pp are assumed to be the same for all combinations, and Tgm is in units of the critical temperature for a 
^Li-^Li mixture. 

function tends to a constant value, F{k oo) = The critical temperature is then given by 

Tgm = 22/3e-i/6^iiexp(-l/A) = 0.825 ;i^exp(-7rV2pF|a|), « > 1. (21) 

Note that this result is quite different from the BCS critical temperature of Eq. ^ . Aside from a constant of the 
order of unity, it contains an extra small factor 1/Vk ~ m/M <C 1. Thus, for a fermionic mixture with a large 
mass ratio of the components the second order contribution significantly reduces the critical temperature compared 
to the prediction of the simple BCS approach. 

In Table 1 we show the critical temperature Tgm following from Eq. (|19p with F{k) given by Eq. ((20l) . for various 
mixtures of fermionic atoms. The critical temperature is given in units of Tqm for a ^Li-^Li mixture, and it is assumed 
that the quantity exp(— 7rS/2pi?|a|) is the same for all mixtures. One clearly sees that replacing one species in 
a mixture by a lighter one increases the critical temperature, whereas the replacement with a heavier one decreases 
Tgm- 

IV. CRITICAL TEMPERATURE. HIGHER ORDER CONTRIBUTIONS 

Let us now consider the contribution of higher order (~ A^) many-body corrections. As can be seen from Eq. (|18p. 
these corrections enter the exponent for the critical temperature being divided by and, therefore, the corresponding 
term is \Q{M / ra) . where Q is a function of the mass ratio. For moderate values of M/m, the function Q is of the 
order of unity, and, therefore, the corresponding corrections can be neglected. However, as we will see later, for 
a large mass ratio, the function Q becomes proportional to M/m^ and the related corrections in the exponent are 
^ kpaM /m. The applicability of the perturbation theory requires {kFofM/m < 1 (sec Eq. ^ below), but the 
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quantity kpaM/m should not necessarily be small for M/m ^1. As a result, the third-order many-body corrections 
proportional to M/m have to be taken into account. 

For calculating the terms of the order of A^^, we have to modify the gap equation ([2]) in order to include the difference 
between particles and quasiparticles, or single-particle excitations [38[. In the considered case of a two-component 
Fermi gas, the quasiparticles are characterized by the effective masses M* and m* and by the Z-factors Zm and Z^. 
The Z-factors are related to the amplitude of creating a quasiparticle by adding an extra particle to the system (see 
[l| for rigorous definitions and details). For a given (effective) interaction Vcs between heavy and light fcrmions, the 
interaction between the corresponding quasiparticles is simply equal to ZMZmVeS- It is important for our approach 
that the ratios M* /M, m* /m, and the constants Zm, Z.,n differ from unity only by a small amount proportional A. 

We can now extend the analysis of the gap equation to higher order terms. As follows from the previous discussions, 
the qualitative form of the gap equation can be written as 



I = V*pZMZmVcS 



Inii 



C 



(22) 



where v*p — 'm*kp /n'^h? and the reduced mass m* is determined by the effective masses M* and m* . The effective 
interaction Vctt = g + SV + SV^'^'' now includes also the third-order many-body contribution SV We are interested 
in the terms of the order of A^ and A"^ In n/Tc on the right-hand side of Eq. ([22|) . where the A"^ In^/Tc contributions 
come from the integration over momenta near the Fermi surface, whereas the momenta far from the Fermi surface 
result in A^ contributions. The term of the order A^ in Eq. ((22|) comes only from the second-order term SV in the 
effective interaction. For a large mass ratio this term contains only ln(M/m) for large mass ratio (see. Eqs. (fT5|) 
andfTT]) and, therefore, can be neglected. The term A^ In^/Tc results from the third-order term SV^^"^ in the effective 
interaction (with VpZ^iZm — > lyp) and from the difference between particles and quasiparticles, {i/pZM Zm — vp) ^ A^, 
multiplied by the first order term g in the effective interaction. As a result, up to terms of the order of A^, we can 
write the linearized renormalized gap equation as 



A(P) 



2Trh^a ml 



-L 



p<Apf 



mr 
dp' 



-ZmZ. 



dp' 



tanh[ei(p')/2Td + tanh[e2(p')/2Tc 



SV{p,n'pp) + 

, n'pp) 



2[a(p') + 6(p')] 

tanhKi(p')/2r,: 



2mr 



A(nV) 



tanh[6(p')/2T, 



2[ei(p')+6(p')] 



A(nV) 



(23) 



As in Eq. (|12p . we introduce an upper cut-off App for the purpose of convergence of integrals at large momenta. Eq. 
([23|l can be solved in the same way as Eq. (fT2|) and wc obtain 



Tc^-e'' y/fiifi2 cxp ■ 

TT 



Tbcs exp >j - ^ 



TcGM exp <{ - — 



I* f 



-^ZuZ„^ - \\ \^vp[W + (5y(3) 




-LZmZ„, - \ \\^vpWi'i) 

, mr- 



(24) 



Note that only the contributions that are linear in M/m for M/m 3> f should be kept in m^/mr, Zm, Z^, and SV^'^K 
The effective masses M*, m* and the constants Zm, Z^n can be obtained from the derivatives of the corresponding 
self-energies S m (^ , V) and S„i (w , p) with respect to the frequency lo and momentum p, evaluated at w = and p = pp: 



ZM{m) — I 1 ^ 




M(?i 



UJ = 0,p=pF, 

M d^M{iU,p) 



pp dp 

m aS]„(ci;,p) 
pp 



dp 



u—G,p—PF) 



UJ—0,P—Pf 



(25) 
(26) 
(27) 



The diagrams for the self energies E^/ and E„j up to the second order in g are shown in Fig. [3] The corresponding 
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Figure 3: The first and second order diagrams to the self-energy T;a{uj,p). Note that we present here only irreducible diagrams, 
and, therefore, the diagram of the second order containing the first order self-energy insertion to the Green function of the 
/3-fermion in the first diagram is omitted. 



analytical expressions read 



^ , , 2TTh'^a 
rrir 

duji dpi 



■6/3(^^2, P2) 



27r {i-Khf 

Gp{uji + W2,P2 +Pi)Ga(-wi -Ka;,-pi + p) - G^"'(wi, pi)G(,"'(-wi, --pi) 



27r {2TThf 

where a = M, (3 ~ m or a ~ m, f3 ~ and the Green functions are given by 

Ga(/3)(w,p) 



(28) 



Ca(/3)(P) +«'5sign[^„(m(p)]' 



<^i%)('^'P) = 2a5 — -7 

with S = +0. The divergent integral in the second-order contribution is renormalized in a standard way by replacing 
the coupling constant g with the scattering amplitude a and subtracting the product G^^\uJi^pi)Ga\~OLJi, —pi) 
of the two Green functions in vacuum {fJ.1.2 = 0), which corresponds to the second order Born contribution to the 
scattering amplitude, from the integrand. After integrating over the frequencies uji and L02 in Eq. (|28p . we obtain 

27r^ ^ f27rfi'aY f dp^ dp, / (1 - /K/3(P2)])/Ka(-Pi + p)]/K/3(P2 + Pi)] 



rrir \ riir J J (27rft)3 (27rft)3 [ w - ^^(-pi + p) - C/3(P2 + Pi) + C/3(P2) - 

/[C/3(P2)](l-/[ea(-Pl+P)])(l-/K/3(P2+Pl)]) ^ fMP2)] 
^ - £.a{-Pl + P) - ?/3(P2 + Pi) +^/3(P2) + i5 p\/2mr - iS 



(29) 



As we discussed above, for a large M/m only the leading contributions that are linear in M/m should be kept, and 
lengthy calculations with the use of Eqs. p5|) - ((29)) give 

Zm = 1, (30) 

1^ -1+5 (^1^2- ^K^) ^' (31) 
Z^ = l-\il + 2ln2)(^y'-i, (32) 

— = ! + -(! + 21n2)M^ — . (33) 
TO 3 V TTrt / m 

The diagrams for third-order contributions to the effective interaction V;ift(p, p') are shown in Fig. 21 where we 
omit the diagrams that can be obtained by inserting the first order self-energy blocks (the first diagram in Fig. [3]) 
into the internal lines of the second-order diagram from Fig. [1] . These self-energy contributions (the first term in Eq. 
([28)) ) simply shift the chemical potentials. It turns out that only diagrams a, b, and c could contain terms linear in 
M/m, whereas the rest of the diagrams are proportional to ln(M/TO). The divergencies at large momenta in diagrams 
d and f can be removed by renormalizing the coupling constant g in the second-order diagram in Fig. [TJ Analytical 
expressions for the diagrams a,b, and c are the following: 

^K^^^(P,P') 



10 



6V^'\p,P') = 



P 



-p. 



+ 



-P ■ • * • ■ -P' 



+ 



-p' -p 



— — « 

J 


» — ^ 


b 







-p' 



p - m — • • - P' p > • y f > p' 



-p ■ • 4 4 ■ -p' -p I » « l i I -p' 



Figure 4: The third-order contributions to the effective interaction between heavy (thick line) and light (thin line) fermions. 
The dashed line corresponds to the coupling constant g. 



dujidpi 

(271)4^3 



Gj\/('^i,P + pi)GM(t^i,p' + Pi) 



duj2dp2 



Gm(w2,P2)G„i(tJi +t^2,Pl +P2), 



6V^^'\p,p')=g' I ^^^ii^G™(c^i,p + pi)G™(coi,p' + pi) / ^5^^Gm(cc>2, P2)Gm(coi + 0.2, Pi + P2), 



where q = p — p'- The integration over frequencies uji and lu2 in the above expressions is straightforward and gives 

e,.(3). „M 3 f dpi /Kl(Pl+q)]-/Kl(Pl)] f dp2 /K2(P2+q)]-/K2(P2)] 



3 1 



(27rft)3 ei(Pi+q)-ei(Pi) J i2nhr 6(P2 + q) - 6(P2) 



where z/j\f = Mpf/2t: h and i^,,,. = mpp /^ir h are the densities of states at the Fermi level for heavy and Ught 
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fermions, respectively. For the other two contributions we obtain: 

/K2(PI+P)I 1 



e2(Pi + p)-6(pi + p') + *'5 e2(Pi + p) 

i-/[6(Pi + p)] 1 



6(pi +p) -6(pi + p') 6(pi+p) + s 



with 



and zero otherwise. A similar expression is obtained for SVc (PjP')? with the replacements ^2 ^ ^1 and Am{s,p) 

Ar,^{s,p). 

The corresponding contributions to the s-wave scattering channel can be obtained by averaging over the directions 
of the momenta p and p': 

<^=/S/f<'(P.p'), . = <..<.,c. 

In the limit of Af/m ^ 1, the leading terms in these contributions are 



,y7>,,^±imrj^) ,35, 

16 \ nh / m 



2 M 



dV^^^^ = 0, (36) 



7-757 l + 4(2 + 31n2)ln2 /apF\2 Af 

dVc' = -g — ( — — , 

18 \ nh / m 



(37) 



where C(a;) is the Riemann zeta-functions (C(3) = 1.202). As a result, the quantity SV^^'f in Eq. (IM)) is 



,51/(3) = ^1^(3) ^ gy{3) ^ gy{3) ^3g^ 



Note that the validity of the perturbation theory requires the quantity (JV^t^) smaller than the coupling constant 
g. This leads to the condition 

{apF/hfM/m < 1. (39) 



Thus, the actual small parameter of the theory in the limit of a large mass ratio M/ni is {pp\a\/h) \J Mjra. 

After substituting Eqs. (f30|) -([33 |) and ([38|l into Eq. (j24|) we find the critical temperature in the limit of M/m ^ 1: 

T, = ^2V3e"V6^ exp {-,r^ - 0.034^^) = Tgm exp{-0.011(p^^|a|/ft)A//m}. (40) 
TT ZM \ Z\a\pF Trn. m / 



Compared to the transition temperature in the GM approach, Ea. (|40p contains an extra exponential factor which, 
in principle, can be large. However, this requires a very high mass ratio M/m. The extra term in the exponent of 
Ea. (|40|) can be written as Q.Ol^J M /m x [pF\a\/h)yJ M /ra and. since the second multiple in this expression is small, 
one should have the mass ratio at least of the order of thousands in order to get a noticeable change of Tc compared 
to the GM result. In this case Eg. ([55)1 shows that the parameter A = 2pp\a\/nh should be very small and, hence, the 
transition temperature itself is vanishingiy low. 

We thus see that for reasonable values of pF\a\/h satisfying Ea. ([39)) . let say p_F|a|/fi ^ 0.1 and M/m < 100, the 
higher order contributions do not really change the GM result for the transition temperature. 

At the same time, our analysis shows that for M/m ^ 1 the conventional weakly interacting regime requires much 
lower values of A = 2pp\a\/TTh than in the case of equal masses and the small parameter of the perturbation theory 
is given by Eq.([39]) In the next section we discuss the physical origin of this parameter. 
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V. 



SMALL PARAMETER OF THE THEORY 



There are several conditions that allow one to develop a perturbation theory for a many-body fermionic system 
on the basis of Hamiltonian ([T|). First of all. this is the condition of the weakly interacting regime, which assumes 
that the amplitude a of the interspecies interaction is much smaller than the mean separation between particles. The 
latter is of the order of h/pp, and we immediately have the inequality 



At the same time, inequality (|4ip allows one to use the binary approach for the interparticle interaction. Then, 
assuming a short-range character of the interatomic potential, one can consider the interaction between particles as 
contact and write the interaction part of the Hamiltonian as g J (ir-0j''(r)-0i(r)'0^(r)?/;2(r). 

In the weakly interacting regime only fermions near the Fermi surface participate in the response of the system to 
external perturbations. Therefore, there is another condition that is needed for constructing the perturbation theory. 
Namely, we have to assume that for both light and heavy fermions the density of states near the Fermi surface is not 
strongly distorted by the interactions. This is certainly the case if both Fermi energies, p|,/2m and p|,/2M, greatly 
exceed the mean- field interaction ng. For M ^ m this condition is equivalent to inequality (|^T|) . In contrast, for 
M ^ m Eq. (|4T|) only guarantees that the Fermi energy of light fermions is p\/2m ^ ng, whereas the condition 
p|,/2M ^ ng leads to the inequality (pi?|a|//i)M/m ^ 1. This mean-field condition, however, is far too strong 
because at the mean-field level, the interaction shifts uniformly all energy states and, hence, results only in the change 
of the chemical potential. Actually, the interaction-induced modification of the density of states is determined by 
the momentum and frequency dependence of the fermionic self-energy (see Eqs. (|25p - (|27p ) and appears in the second 
order of the perturbative expansion in g (the second diagram in Fig. [3]). The corresponding contribution describes 
the process in which a heavy fermion pushes a light one out of the Fermi sphere and then, interacting once more 
with this light fermion, puts it back to the initial state. Due to the Pauli principle, the momenta of both light and 
heavy fermions in the intermediate state should be larger than the Fermi momentum. As a result, for the initial 
heavy-fermion state close to the Fermi surface, the most important intermediate states will be those with momenta 
close to the Fermi momentum. Therefore, the resulting contribution should be proportional to the product of the 
densities of states of heavy and light fermions at the Fermi surface, and the relative change of both densities of states 
is controlled by the parameter g^vuVm ~ {ppa/Kf'M /m. Thus, this parameter should be small, i.e. we arrive at 



A complementary physical argument on support of this small parameter comes from the consideration of the effective 
interaction between a light and a heavy fermion in the medium. For example, the process described by the diagram in 
Fig. 4a can be viewed in the following way. Incoming heavy and light fermions interact with fermions inside the filled 
Fermi spheres and transfer them to the states above the Fermi surfaces. Then the transferred heavy and light fermions 
interact with each other and return to their initial states. The important point is that the intermediate state of this 
process contains excitations (particle-hole pairs) near the Fermi surface of the filled Fermi sphere of heavy fermions. 
Therefore, the corresponding contribution to the effective interaction is geS ~ g'^VmVM, where the densities of states 
of heavy and light fermions near the Fermi surface are vm = Mpp /2TT'^hP and = mpp /'iir'^h? , respectively. This 
leads to geS ~ g{pFa/h)'^M/m. Comparing it with the direct interaction g and requiring the inequality |geff| ^ \g\ 
which allows one to use a perturbation theory, we again obtain a small parameter of the theory {ppa/h)'^M/m <C 1. 

Let us now understand in which physical quantities the parameter (|39p enters directly. In the limit of M/m ^ 1 
heavy fermions occupy the energy interval p']p/2M which is much narrower than the energy interval p\/2m occupied 
by light fermions. However, the heavy-fermion density of states is much larger: vj^i — Adpp / (2t^'^1t?) 3> Vm = 
mpF / (2Tr^h^). The high density of states of heavy fermions manifests itself in any quantity characterized by processes 
where, for the heavy fermions, only the states near the Fermi surface are important. This is the case for the effective 
masses of atoms, critical temperature, and (see the next section) for the zero-temperature order parameter Aq. If, 
however, all states of the heavy fermions are important, then the peak at i? ~ p\/2M in the energy distribution 
of heavy fermions is integrated out, and the result does not contain the parameter (I39p . This is exactly what is 
happening in the calculation of the second order correction to the energy of the system, which involves the sum over 
all energy states. 

The second order contribution to the energy is shown diagrammatically in Fig. [5l and the corresponding analytical 
expression reads 



PF\a\/h < 1. 



(41) 



Eq. jSSl): 



{pFa/hfM/m <C 1. 




(42) 
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t 




Figure 5: The second order contribution to the energy of the system. 



where 



duji dpi 



Ga (w + wi , q + pi )Gq (tJi , pi ) 

/Ka(Pi+q)]-/Ka(Pi)] 



{2TThf u - (Ca(pi + q) - Ca(Pi)) + i(5(sign[$a (pi + q)] - sign[^„(pi)]) 



is the polarization operator (the bubble in the diagrammatic language) for a-fermions [a — M or to). As it can be 
seen, the integration over oj in Eq. (|42p results in an integral that diverges at large q. This divergence, however, can 
be eliminated by subtracting the second-order Born contribution to the interparticle scattering amplitude multiplied 
by the densities of fermions. This corresponds to the renormalization of the coupling constant g in the first order 
(mean-field) contribution to the energy E^^'^ = g-n? (see [ij for more details). The resulting expression then coincides 
with equation (6.12) in 

After using the spectral representation for the polarization operator Ha{'^,p)'- 



Jo 



1 



1 



with the hmction Aa{s,p) from Eq. (j34p . equation 



uj — s + iS u} + s — iS 
can be rewritten in the form 



i5(2) 



dq 



'ds, r d,/M{s^,q)A^{s,,q) 

Jo Si + S2 



2mr 



(43) 



where the second term in the brackets corresponds to the renormalization. In the limit of M/m ^ 1, as follows 
from Eq. (|34p . typical values of si are much smaller than typical values of S2. We therefore can neglect si in the 
denominator of the first term in Eq. (j43p and replace to^ by to. This gives 



2-^ 



= 2' 



dq 



2^ 



{2TThf 

dq 

{2-Khf 
dq 

{2TThf 

dq 



dsiAM{si,q) 



ds- 



Am{s2,q) 2m 



S2 



dsiAM{si,q) 



1 



n™(o,(7) 



9{q - 2pf) 



gn 



{2Trh)3 
2 9(81n2- 9) apF 



0{q - 2pf) 



ApF 



Up- 



2m n 
q^ 

{2pF- 



12pl 



e{2p 



F 



q) 

-q) 



n™(o,(7) 



1 + ^ 



2to 



In 



2pi 



2to 



140 TTh 

As we see, the final result does not depend on A/, as it was anticipated above. 



VI. ORDER PARAMETER AND SINGLE-PARTICLE EXCITATIONS 



We now calculate the order parameter and its temperature dependence. At zero temperature no quasiparticles are 
present since both quasiparticle energies are positive {E± > 0). Then, confining ourselves to second order terms in 
the renormalized gap equation reads: 
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Ao(p) = - 



2Trh'^a 
rrir 
dp' 

{2TThy 



dp' 

(2^ 



-SV{p,p') 



1 



E+{p') + E-ip') 
1 



2771,. 
„/2 



Ao(p') 



i?+(p') + i?-(p') 



7tAo(p'). 



(44) 



Strictly speaking, the many-body contribution to the interparticle interaction SV is affected by the superfluid pairing 
and, therefore, does not coincide with that of Eq. (jlip . However, at zero temperature the difference is proportional 
to Aq /fj,i and, hence, is exponentially small. Therefore, we can use Eq. ((TT|) for SV in Eq. ([44]). 

The arguments used above for obtaining Eq. from Eq. can also be applied here, but the large logarithm 
ln{iJ,/Tc) should be replaced by ln(/i/Ao(pF)). 

For the value of the order parameter at the Fermi surface, Ao(pf)j we obtain 



Ao(pf) 



7/3 



exp 



1 + ln4 



[/(«) + /(«-!) -1] 



4777,,. 



exp I 



(45) 



Comparing Eq. ^TE\\ with Eq. ps)) we obtain a relation between Ao(pf) and Tc 



TT K 



1/2 + 



T7iAo(pf). 



(46) 



We should emphasize that relation (|45p between the order parameter and the critical temperature remains valid after 
taking into account higher order terms in the gap equation, which is necessary for a large mass ratio (see the previous 
section for the discussion of the critical temperature). The generalization of Eq. (j44p in order to include the higher 
order terms repeats the derivation of Eq. (j23p and the resulting equation reads: 



A(P) 



27rfi,^a ; 



rrir 777, 

dp' 

(277/7)3 



-ZAfZ„ 



dp' 



(277/7)3 



1 



2mr 



E+ip') + E^ip') 



A(nV) 



p<ApF 



6V{p, n'pp) + 



E+{p') + E^{p') 



7TA(nV)- 



This equation can be solved in a way similar to that of solving Eq. (j23p , and the solution is 



(47) 



Ao(pf) = TTe ^ 



A. 

Arrir 



exp ' 



-ZmZ^X - i/F (sV + (51/(3)^ 



(48) 



Comparing Eq. p4|) with Eq. (|48)) . we immediately obtain Eq. (|46l) . 

We now analyze the order parameter in the two limiting cases: M = m and M ^ m. In the case of equal masses 
we have 777^ = 777/2 and recover the usual expression [30| for the order parameter from Eq. (|45p : 



Ao(pf) = ( 7 j TTTT^xp 



2M 



2pF\a\ 



0.489 Tc 



A/ = 777. 



In the case of M 3> 777, i.e. k 3> 1, we have 777^ ~ 777 and, including higher-order contributions, from Eqs. (|46p and 
21\i we obtain 



Ao(pf) = 2«/3e-i3/6 PLexp f- 
2v M777 V 



2pF\a\ 



0.011 



|a|pF M 



0.882 



> T^; M > 777. 



(49) 



Note that in this limit the order parameter at the Fermi surface is much larger than the critical temperature. 

In order to analyze the behavior of the order parameter A for temperatures close to the critical temperature, 
(Tc — r) <C Tc, we have to expand the gap equation ^ in powers of A/Tc <C 1 and keep the cubic term. The result 
can be written as 



Ak 7C(3) M(pF) y 
r (1 + k)^ 87r2 { J 




Figure 6: Two branches E±{p) of single-particle excitations for M/m = 10. 



From this equation we obtain 



In the limit of equal masses, k = 1, we reproduce the well-known result for the temperature dependence of the order 
parameter. In the opposite limit of a large mass ratio, k ^ 1, we find 



Note that as well as Eq. the obtained equation (|50p contains a large factor ^jM/m. 

At any temperature the energies of single-particle excitations are given by (see Eq. Q and Fig. [6]) 




i?±(p)=± £— ^ +W ^ +A2(p^), (52) 



where m_ = Mm/{M — m). Equation (|52p reveals a peculiar feature of fermionic mixtures with unequal masses of 
components. For equal masses, the minimum of -E'±(p) (i.e. the gap) is at the Fermi surface and equals A(pi?). The 
situation for unequal masses is different: the single-particle excitation energies E±{p) reach their minimum values 
E±rnin = A{pf)VI - s2 = 2A{pf) / (k'^^'^ + K''^^'^) at momenta p^f = p|,zp2mKi/2(K-l)A(pi^)/(K+l). For M/m > 1 
the corresponding gap is much smaller than A{pp): 

E±mUM » m) « 2A(pf )«;-i/2 ^ 2J^A{pp) « A{pp). (53) 



It is interesting to note that the presence of a small factor -y/ m/M in Eq. ([53|) restores the intuitive picture that 
the gap in the single-particle spectrum and the critical temperature are of the same order of magnitude even in the 
limit of a large mass ratio. We point out, however, that in this limit the order parameter on the Fermi surface, being 
much larger than the critical temperature, is not equal to the gap in the single-particle spectrum. This gap is of 
the order of the critical temperature, and the low-energy single-particle excitations correspond to momenta different 
from the Fermi momentum pp. Owing to the former circumstance, one does not expect any dramatic changes in 
thermodynamic properties of the system with increasing the mass ratio M/m to a large value. 
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Figure 7: The contribution of three-body processes between one light and two heavy fermions described by the connected 
three-body vertex ri"^' , to the Gorkov-Melikh-Barkhudarov corrections. 



VII. THREE-BODY RESONANCES 



Let us now discuss the influence of the three-body physics on the results of the previous sections. The diagram for 
the Gorkov-Melikh-Barkhudarov corrections in Fig. [1] corresponds to collisions between three particles: two from a 
Cooper pair and one from the filled Fermi sea. Thus, it is a three-body process. During this process, however, the 
three particles undergo two successive two-body collisions and never appear simultaneously within the range of the 
interatomic interaction. The corresponding three-body wave function vanishes when the hyperspherical radius (see 
the definition before Eq. ([51]) ) is tending to zero. 

In a dilute two-component Fermi gas, real three-body collisions during which the three colliding particles simul- 
taneously approach each other, are rare. An additional smallness compared to a Bose gas is provided by the Pauli 
principle. Two of the three colliding particles are identical fermions and, therefore, the wave function of their relative 
motion should strongly decrease at small separations. As a result, the contribution of such collisions to the effective 
pairing interaction is small and can be neglected. However, for the case of a large mass ratio M/m the situation is 
more subtle. If M/m > 13.6, two heavy and one light fermions can form three-body bound states [23,[3^li3|- The 
most interesting case corresponds to the presence of a weakly bound trimer state because this results in a resonance 
3-body scattering at low energy. It is not clear that these resonances should be taken into account when calculating 
the pairing energy since there may be nontrivial issues of wave function statistics involved; nevertheless we shall 
estimate their possible contribution and leave such issues for further study [s^ . 

To analyze the effect of three-body bound states we note that the contribution to the Gorkov-Melikh-Barkhudarov 
corrections in Fig. [T] is part of a more general contribution involving the connected three-body vertex function 
ri'^^(see Fig. [7]). This is a consequence of a general relation between two- and three-particle Green functions. The 

quantity Fj; ({Pilin ' Wlout) '^^^'^ Pi = i^ii pOj * = Ij 2, 3, describes the scattering of two heavy and one light particle 
from the initial state with incoming energy- momenta pi into the final state with outgoing energy-momenta p'^. For 
LUi = pj /2m (the mass-shell condition) the vertex function coincides with the T-matrix. By definition, the connected 
vertex function Fc does not include three-body processes in which only two out of the three particles collide (in 
our case, a light fermion collides with only one heavy fermion) and, therefore, Fc is represented only by connected 
diagrams. The general three-body vertex function F'^'^)(see Eq. [56)) contains all diagrams including disconnected ones. 
Those describe processes in which only two out of the three particles interact with each other. Fig. [5] shows the 

(3) 

simplest contribution to Fc that is second order and results in the Gorkov-Melikh-Barkhudarov corrections shown 
in Fig. n 

We consider the case where the size of a three-body bound state is much larger than |a|, but much smaller than 
the average distance between particles in the gas. Accordingly, the binding energy is much larger than typical kinetic 
energies of particles, the Fermi energies Hi. In this case, the influence of other particles of the gas can be neglected 
and the properties of the bound state can be found by solving the three-body Schrodinger equation. Introducing the 
hyperspherical radius p = •\/x^~+"y^ in the 6-dimensional space (x,y), where x-\/ (2M -I- m) /4m = ri — (r2 -I- r3)/2 
is the distance between the light fermion and the center of mass of two heavy ones separated from each other by a 
distance y = r2 — ra, the normalized wave function of a shallow bound state with the binding energy Ei, = —1t?'/2MIP' 
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Figure 8: The lowest order contribution to the connected three-body vertex. 



and the size b ^ \a\, has the form [41| 



g , a p 



Here $i and $2 are the functions of hyper angles fl, and we do not give exphcit expressions for these functions 
because of their complexity. For p > b, the wave function decays exponentially. Note that the normalization of the 
wave function ([54)) is determined by distances p ~ \a\. 

Most conveniently the contribution of the bound state to the vertex function can be found using the three-body 
Green function G{{pi} , {p-} , w): 



uj-H + iO 



where the Hamiltonian H has the form H = Hq + V with 



and 

V = g5{Ti - V2) + g5{vi - rg) 

in the coordinate representation (index 1 corresponds to the light fcrmion and indices 2 and 3 to the heavy ones). 
The Green function satisfies the equation 

7?G(3)({pj,{pa,^)= n ^(p^-p^' 

i=l,2,3 

which is equivalent to the integral equation 

G(3) =G^=^)+G^')i/G(3), (55) 

with Gq'^'' = [oj — p\l2m — (jr^ + P3) /2A/ + iO]^^ (5(pi — p^) being the Green function for free particles. This 
equation can be rewritten in the form 

G^^^ ^Gf^ +Gf^T'^^^Gf\ (56) 

where we introduce the vertex function T^^\ This function describes all scattering processes involving three particles, 
both connected (described by Fc^^) and disconnected ones (not included in Fc'^^). The vertex function F^'^^ obeys the 
Lipmann-Schwinger equation 

f(^^) = 1/ + T>g[,^^f(3) (57) 
and, as it can be seen from Eqs. ((55 |) -([57 |) . is related to the Green function G*-^-* as 

f(3) = y-f i/G^^^y. (58) 
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It is convenient to use the spectral decomposition of the Green function. In the center-of-mass reference frame, where 
Si Pi — X)i Pi ~ 0, this decomposition reads: 

G(='^({pJ,{p:},^) = E*»({P^}) -H*"({Pa)+ / d\¥^>i{p,}) l^*(+)({p^}), (59) 

OJ — £jn + lO J OJ — hx + lO 

where the summation is performed over a complete set ^'^n, ^a^"* | of eigenfunctions of the three-body Hamiltonian H 

with eigenenergies E\, respectively. The eigenfunctions correspond to bound states with energies En < 0, and 

the eigenfunctions to scattering states of three particles with energies Ex > 0. Their asymptotic behavior at large 
interparticle distances contains incoming plane waves with momenta specified by the index A, and outgoing (therefore, 
index +) scattered waves. These eigenfunctions vanish for small hyperspherical radius p and, in particular, they 
describe the Gorkov-Melik-Barkhudarov corrections discussed above. On the contrary, the bound state eigenfunctions 
are nonzero for small p and decay exponentially for p ^ oo. 

Eqs. (|58p and ((59|) together give the decomposition of the vertex function F'^^ in terms of the solutions of the 

three-body Schrodinger equation. Obviously, the bound states contribute only to the connected vertex function Tc ■ 
In particular, the contribution of the bound state ^„ is: 

5„r(3) ({pj , {p^} , u) = 5„r(3) ({pj , {pa , ^) = p*„({p,})l * L^v-^niW,}). m 

By using the Schrodinger equation 

(iJo + V)^n = En^n, 

Eq. (I60p can be rewritten in the form 

^nr(3) ({p.} , {pa , = ^:({pa) ' ^""p-""^- -mwm 

U) — En + «0 

where i^o({Pz}) = Pi/2m + {pl+pf} /2M. 

Now we can estimate the contribution of the weakly bound state of one light and two heavy fcrmions to the effective 
interparticle interaction Vcs between light and heavy fermions with opposite momenta on the Fermi surface. The 
analytical expression corresponding to the diagram in Fig. [7] is 

^K^ff = / ^ / 7;5%r(3) (p, q, -p; p', -p', q; u) . (62) 



The contribution of the bound state is then obtained by substituting Eq. (1611) into Eq. (|62p and integrating out the 
frequency lo: 

J [2Trhy [q^ - pp)/2M + pi + 2p2 - En - q'^/[2{2M + m)\ 

where the last term in the denominator corresponds to the motion of the center of mass. In the considered case, the 
binding energy En is the largest energy scale {En S> pi,p2)i and SnYcff in Eq. (|63p can be estimated as 

SnV,e ^ {pF/hfEb\^n{0,0,0)\\ (64) 

where the factor {pp/Kf' results from the integration over dq and we used the condition \a\pp/h ^ 1 to set all 
momenta in the wave function of the bound state to zero. After using the wave function from Eq. (|54p . we obtain: 



3/4 



2 1,2 



and therefore 



*„(0, 0, 0) ~ {pp/h) ( — 1 {pp/h)a'b 



vpSnV^s ^ MPF/hfEba%'' — ^{pp \a\ /h)\ppb/hy 

1 ™ ; " m 
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Figure 9: Regimes of superfluid pairing for M/m 2> 1. In the intermediate regime the conventional perturbation theory is not 
apphcable (see text). 



This result has to be compared with the GM contribution vpSV ^ [pp |a| /Kf' ln(M/m), and with the contribution of 
third-ordrer terms vp5V^^^ ^ {pF\a\/hYA4/m. Under the condition {ppb/h) < 1 corresponding to a not too shallow 
bound state, we find that the contribution of three-body resonances is small compared to both of them: 



=^ - {PF \a\ h) (pFb h) < 1, 

I'pSV ln(i\7/m) 



and 



ipF\a\/h){pFb/h)^^m/M <^1. 



We thus see that three-body resonances are rather narrow, and their contribution to the effective interaction can 
be omitted. So, the results obtained in the previous sections for the critical temperature, effective masses, order 
parameter, and elementary excitations remain unchanged. 



VIII. CONCLUDING REMARKS 



We now give an outlook on the physics of attractively interacting mixtures of heavy and light fermionic atoms 
in view of the results obtained in this paper. We have developed a perturbation theory in the BCS limit for the 
heavy-light superfluid pairing along the lines proposed by Gorkov and Melik-Barkhudarov [s^l and found that for 
M/m ^ 1 one has to take into account both the second-order and third-order contributions. The result for the critical 
temperature and order parameter is then quite different from the outcome of the simple BCS approach. Moreover, 
the small parameter of the theory is given by Eq. ([5^ and reads: {pp\a\/h) ^ 1. As we explained in Section V, 
this can be seen from the second-order correction to the fermionic self-energy, which is controlled by the parameter 
g^vj^iVm ~ [ppa/hY'M jm. Therefore, in a mixture of heavy and light fermions the conventional perturbation theory 
for the weakly interacting regime requires much smaller pp (densities) and/or \a\ than in the case of M ^ m, where 
the small parameter is {pp\a\/h) <C 1. 

Let us now discuss the cases of M = m and M ^ m regarding the regimes of superfluid pairing. For M ~ m we 
have the strongly interacting regime for {pp\a\/h > 1, and the BCS limit for {pp\a\/H ^ 1. In the former case the 
perturbation theory is not applicable and the results are obtained either by Monte Carlo methods or by adjusting 
the mean- field theory to this regime (see Q for review). For M ^ m the situation is different. As we found, the 
conventional perturbation theory works well under the condition {pp\a\/h) <^ yJmjM (see Fig. |9]). For {j)p\a\lh > 1 
we have the strongly interacting regime where the perturbation theory does not work at all. However, we now have 
a range of densities and scattering lengths, where ^ m/M ^ {pp\a\/h) <C 1. In this intermediate regime one can still 
use Hamiltonian ([T]) and try to develop a perturbative approach, since the scattering amplitude is much smaller than 
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the mean interparticle separation. On the other hand, the conventional perturbation theory does not work for the 
reasons explained in Section V. In order to construct a reliable theory one should at least renormalize the interaction 
between heavy and light fermions by making an exact resummation of diagrams containing loops of heavy and light 
fermions. We then expect a substantial renormalization of the properties of the superfluid phase. 

We thus see that our findings pave a way to revealing novel types of superfluid pairing in mixtures of attractively 
interacting ultracold fermionic atoms with very different masses. An appropriate candidate is a gaseous mixture 
of ^^^Yb with ^Li, and one should work out possibilities for tuning the Li-Yb interaction in this system. Another 
candidate is a two-species system of fermionic atoms in an optical lattice with a small filling factor. The difference in 
the hopping amplitudes of the species can be made rather large, which corresponds to a large ratio of the "heavy" to 
"light" effective mass. For example, in the case of ^Li-'*°K mixture one can increase the mass ratio by a factor of 20 
in a lattice with period of 250 nm and the tunneling rates ~ 10'^ s~^ and ~ 10^ s~^ for K and Li, respectively. 
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